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The counting statistics give insight into the properties of quantum states of light and other 
quantum states of matter such as ultracold atoms or electrons. The theoretical description of photon 
counting was derived in the 1960s and was extended to massive particles more recently. Typically, 
the interaction between each particle and the detector is assumed to be limited to short time 
intervals, and the probability of counting particles in one interval is independent of the measurements 
in previous intervals. There has been some effort to describe particle counting as a continuous 
measurement, where the detector and the field to be counted interact continuously. However, no 
general formula applicable to any time and space dependent field has been derived so far. In our 
work, we derive a fully time and space dependent description of the counting process for linear 
quantum many-body systems, taking into account the back-action of the detector on the field. We 
apply our formalism to an expanding Bose-Einstein condensate of ultracold atoms, and show that 
it describes the process correctly, whereas the standard approach gives unphysical results in some 
limits. The example illustrates that in certain situations, the back-action of the detector cannot be 
neglected and has to be included in the description. 



I. INTRODUCTION 

The counting statistics of photons have been used to 
identify the quantum state of light since the beginnings 
of quantum optics. Despite early attempts [lH3[ to use 
a semiclassical approach, the theoretical description of 
the counting process requires a full quantum mechanical 
description of the electromagnetical field and the photo- 
counter. This task was first achieved by Glauber [U-0], 
for an analogous approach see ||, and will be referred 
to as the quantum counting formalism in the following. 
The formalism uses a full quantum mechanical descrip- 
tion of the interaction between the incoming light and 
the atoms in the detector. However, consecutive detec- 
tion events are treated independently, such that the effect 
of absorption of photons at the detector is neglected. 

Even though the quantum counting formalism is being 
applied in a wide range of situations, its limits have been 
pointed out and discussed in the last decades. For ex- 
ample, for long detection times, it may lead to negative 
probabilities or expectation values exceeding the total 
number of photons. Such limitations were first pointed 
out by Mollow [9j and several authors derived a formula 
taking into account the back-action of the detector on the 
field by considering the evolution of the system composed 
of the detector and the field [gMTlj . In 1981, Srinivas and 
Davies [12[ provide a systematic description of photon 
counting as a continuous measurement process, where 
the detector continuously interacts with the light field. 
They derive their formalism based on a one-count oper- 
ator and a no-count operator that determine the time 
evolution of the combined system of the detector and 
the light field. Whereas these operators are postulated 
in the model by Srinivas and Davies (SD model), Imoto 
and coworkers derived a microscopic theory of the con- 
tinuous measurement of the photon number (l3| . They 



derive the no-count and one-count operators taking into 
account the interaction of the photons with the detector 
using the Jaynes-Cumming Hamiltonian for the field of 
a two-level atom. Photon counting including the back 
action of the detector has also been treated in 0, [l5[ . 

Contrary to the quantum counting formalism, the SD 
model predicts that a substraction of one photon from a 
thermal state can increase the expectation value of the 
number of photons [la ]. Recently, this has been con- 
firmed in an experiment where the action of the annihi- 
lation operator on different states of light has been im- 
plemented and measured [lT} • As discussed in [l8| this 
is an experimental confirmation of the SD model. 

The previous discussion suggests, that the SD model is 
a valid description for the photon counting process based 
on continuous measurements. However, a closed formula 
is derived only for the case of a single-mode field and the 
application of the formula to real experimental situation 
is therefore limited 0, [2(|. In 1987, Chmara derived a 
general formula for the photon counting distribution for 
a multimode field [2lj by applying the photon counting 
approach by SD to an open system. While the formula is 
in principle applicable to a wide class of systems, to our 
knowledge, no practical case where the time dependent 
intensities have been calculated has been reported. 

In our work, we study the role of the back-action of 
the detector when counting massive particles, however, 
the results are valid also for photons. Previously we have 
shown that the particle counting statistics can be used as 
a detection method for many-body systems of ultracold 
atoms [22l - |24j . In this work we derive a general formula 
that describes the counting process including the back- 
action of the detector on the field. In contrast to the ear- 
lier descriptions of the counting process stated above, we 
give a fully time- and space dependent description of the 
process. We illustrate the importance of the back-action 
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of the detector on the field by applying our formula to 
the detection of an expanding Bose Einstein condensate 
(BEC) by particle counting. We compare our formalism 
to the quantum counting formalism and to an approxi- 
mate solution obtained by the Born approximation. We 
show that the approximate solutions, although more ac- 
curate than the quantum counting formalism, still fail to 
describe the counting process correctly. 

The paper is organized as follows: In section ITll we re- 
view the quantum counting formalism and the SD model. 
In section IIIII we consider a fully time- and space de- 
pendent description of the counting process and derive a 
formula for the counting distribution which is extended 
for two detectors in Sect. IIV1 We illustrate the space 
and time dependent counting process by considering the 
counting statistics of a freely expanding Bose Einstein 
condensate in Sect. [V] We compare our time- and space- 
dependent counting formalism to the quantum counting 
formalism and show that it some limits, the last one leads 
to a divergent intensity at the detector. We also analyze 
the effect of the absorptive part in the modified field op- 
erator by comparing the exact solution to other approx- 
imative methods. We summarize our results in Sect lVIl 



II. PHOTON COUNTING: STANDARD 
APPROACH AND CONTINUOUS 
MEASUREMENT APPROACH 

In this section, we review the two main approaches 
to photon counting: The quantum counting formalism 
derived in the works of Glauber 0-0] and the continuous 
measurement approach derived by Srinivas and Davies 
12]. 



A. The Quantum Counting Formalism 

The derivation of the quantum counting formalism is 
based on the description of the quantum mechanical in- 
teraction of the photons with the atoms in the detec- 
tor. The approach uses perturbation theory to describe 
the interaction for short intervals of time. The counting 
distribution for the full detection time r is obtained by 
dividing it into small subintervals At and treating the 
measurement in the full interval as a number of succes- 
sive independent measurements in each interval. The ap- 
proach thus describes a sequence of measurements, where 
the field evolves as in the absence of the detector. The 
method is based on the assumption, that the detection 
in one subinterval is independent of the detection in the 
previous subintervals. For the case of a light beam falling 
on a photo detector, it is argued [H, P- 723] that each 
element of the optical field interacts with the detector 
only briefly. The response time of the detector is short 
and the energy of the electron state is well defined after 
an interaction time At. For such a system where the 
unabsorbed photons escape, there is no need to consider 



the measurement back-action. The resulting equation for 
the counting distribution reads 

P ( m ,t,r) = (T : ^ :). (1) 

ml 

Here T and :: stand for time and normal ordering, respec- 
tively. The intensity operator I(t) is defined in terms of 
the positive and negative frequency parts E ± (r, t) of the 
field operators by X(t) = e J r dt' f Q E~ (r, t)E + (r, t)dr, 
where the spatial integral runs over the detector area f2 
and e denotes the efficiency of the detector. The normal 
ordering reflects the fact that the photons are annihilated 
at the detector. For a single mode field, the formula reads 

p(m,T) = r)(err(eM-^)r- m (n\p\n). (2) 
where p is the density matrix of the measured system. 

B. The Back-action of the Detector on the Field 

In [l2j Srinivas and Davies developed an approach to 
photon counting based on continuous measurements over 
an extended period of time. Their work was motivated by 
the fact that the quantum counting formula Eq. ([1]) ex- 
hibits inconsistencies in the limit of long detection times, 
such as negative probabilities or mean particle numbers 
that exceed the total number of particles. For a single 
mode field the counting probability distribution is given 
by 

p(m,t,T)=J2 ryi~e-n m (zM-cT)) n - m {n\p\n). 

n—m ^ ' 

This expression is formally equivalent to the quantum 
counting formula Eq. ([2]) with the term er substituted by 
1 — e~ er . For er 1, the two formulas are equivalent but 
the SD counting distribution in Eq. §3§ docs not suffer 
from the inconsistencies as the quantum counting formal- 
ism. The mean number of photons that are registered is 
bounded for r —¥ oo and no negative probabilities occur. 
However, the applicability of the formalism derived by 
Srinivas and Davies is limited to a one-mode field and 
therefore fails to describe most experimental situations. 

Summarizing, the quantum counting formalism Eq.Q 
becomes meaningless if er 3> 1, when it can result in neg- 
ative probabilities or unlimited mean number of counted 
photons as r tends to infinity. Most experimentally rel- 
evant situations of photon counting are not in this limit 
and thus are well described by the quantum counting 
formalism. In our work we are interested in the count- 
ing statistics of many-body systems of cold atoms where 
one can control the space and time dependence. In the 
following section, we derive a formula for the counting 
distribution that takes into account the back-action of 
the detector on the field. It generalizes the formalism 
developed by Srinivas and Davies to time- and space de- 
pendent fields. 
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III. PARTICLE COUNTING OF TIME AND 
SPACE DEPENDENT FIELDS 

The time- and space-dependent counting process can 
be described by a master equation that describes the 
interaction between the detector and the detected field 
^(r, t). In typical experimental situations, the interac- 
tion between the detector and the field is restricted to 
a given volume, such as the surface of a photodetector 
or a microchannel plate. We define the function f2(r) to 
describe the spatial configuration of the detector. The 
master equation then reads 

p{t)=e J drQ(r)^(r,t)p¥(r,t) 

-| J drn(r)[*(r,t)t*(r,t)p-p*(r,t)t*(r,t)], (4) 

where p(r, t) denotes the density matrix of the system. 
The first term on the right-hand side of eq. (j4j corre- 
sponds to the number of quantum jumps in the detector 
volume whereas the second term represents the damping 
of the field due to the absorption at the detector. 
In order to solve the master equation Q we first perform 
the transformation pit) = F(t)p{t)F* (t) and define 

*(r,t) = F- 1 (t)*(r,t)F(t), (5) 



where the operator F(t) is defined as 

F(t) = Te~ e/2 ti rf*'/'ir'0(r')* t (r',t')*(r',t') > (g) 

Here, the term T on the left side of the opera- 
tor denotes time ordering, whereas it denotes op- 
posite time ordering on the right side of the op- 
erator. We use the relation e ,A Be~ lA = B + 

2 

7L4, B] + -^r[A, [A,B]] + ... and the commutation re- 
lations for linear fields, [*(r, i), ¥ (r', t)] = <5(r,r'), 
[*(r,t),*t(r',t')] = G (t,r%t'), [¥(r,i),*(r',f)] = 0, 
where G (r, r'|t, t') = (r|e~ l(t_t ' )/riHo |r') is the propaga- 
tor for the time evolution for the Hamiltonian Hq that 
describes the evolution of the field without detection. 

The evolution of the rotated density matrix fulfills the 
equation 

~p{v,t) = e J dr'n(r')V(r,t)p¥(r,t). (7) 

A master equation of this kind, for non-rotated ^ and p 
instead of & and p leads to the quantum counting for- 
malism described by eq. ([1]). Eq. ([7]) can be solved using 
perturbation theory (note that p(0) = p(0)), such that 



p(t)=p(0)+ef*dt'fdr'n(r>)V(r\t')p(0)¥(x',t') (8) 
+e 2 f* dt' f*' dt" J dr' J dr"0(r')0(r")#(a/, t')V(x" ,t")p(0)& [x" , t")¥(x', t') + ... 

I 



We use the cyclic properties of the trace to calculate the 
probability p(m, t) of finding m particles within the de- 
tector opening time r and from the mth order term in 
the expansion in eq. |H]we get 

P (m,T) = (e m J dt' J dt".. J dr' J <Zr"Q(r')ft(r")--( 9 ) 

¥(x', t')¥{x", t")..F^(T)F{T)..y(x", t")$(x", t')). 

We rewrite this expression as a normal ordered expres- 
sion with respect to the modified operators ^, which is 
also normal ordered with respect to the operators $ , as 
they are related by a linear transformation. Taking into 
account the normalization of the counting distribution, 
we obtain 

pKrJHT:*^:), (10) 
m! 

where the intensity at the detector is given by 

I(r) = e^ dt J dvVL{r)¥{r,t)^{Y,t). (11) 



This equation is formally equivalent to the quantum 
counting formula, Eq. ([1]). However, whereas for the 
quantum counting formalism the intensity of particles 
registered at the detector is determined by the square of 
the field operator, in Eq. (f^5|) the intensity is calculated 
using a modified field operator ^(r, t), which includes the 
absorption at the detector. In the following, we analyze 
these modified field operators. 

We rewrite eq. ([5]) by dividing the detection time into 
small sub- intervals At = t/N. The time integration in 
Eq. ([5]) can be written as a sum, such that F(t) = 
U t Fi(fi), with Fi(U) = e -iAt/dr'n(r')*V.**Wr',iO, 
and we get 

*(r,t) = Fi 1 (t 1 )...F N 1 (t N )y(r,t)F N {t N )...F 1 (t 1 ).(12) 

We evaluate the expressions by using the commutation 
relations stated above. We start with the inner term, 

F^(t N )^(r,t)F N (t N ) = e- eAtn W*(r,^). (13) 
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The second term thus reads 

e- eAtn ^F- 1 _ 1 (t N -xMr,t N )F N _ 1 (t N . 1 ) = 



The successive terms are calculated analogously, and in 
the limit of infinitesimal small time intervals we get 



*(r,t) 



dr / G(r,r , ,t,t )*(r',to) J 



where 



G(r,r',Mo) = (r|e- 



-i(t—t a )/h{H —ithtt)\r 



r'). 



(15) 



(16) 



We have thus obtained an expression for the modified 
field operators ^(r, t) which differs from the definition of 
the standard field operator as it includes the propaga- 
tion in an imaginary potential created by the detector. 
Together with the counting formula Eq. (|10[) , this allows 
us, in principle, to calculate the counting distribution for 
time dependent systems with arbitrary detector geome- 
tries. However, solving Eq. (|15[) is in general a highly 
non-trivial task. In Sect. |V]we solve the equation for the 
detection of an expanding BEC. 

It is interesting to point out that in many experimental 
situations, the detection process is fast compared to the 
time evolution of the system. In this case, we can neglect 
the part corresponding to the Hamiltonian H$ in eq. (|15j) 
and get 



*(r,t) = / ^'(rle-^-^V)^!-',^) 

= e- eA <*- to >¥(r,to). (17) 



The intensity ([25]) thus reads 

I(r) - e [ dr(l - e -^ T )¥(r, i„)tf (r, to), (18) 



which is a generalization of the formula Eq. ([3]) consider- 
ing finite detector volumes. For er <C 1, Eq. (fl~8| reduces 
to the quantum counting formula Eq. ([1]) for time inde- 
pendent systems. 



IV. DETECTION WITH TWO DETECTORS 

Our formalism is easily extended to calculate the joint 
counting probability for the detection at two detectors. 
The master equation that describes counting with two 
detectors reads 

p(t) = ei J drill (r)¥(r,t)ptft( r)t ) 

+e 2 J drSl 2 (r)y(r,t)p&(r,t) 

~ J drOi(r)(*(r,i)t$(r,t)p-p*(r,t)t*(r,i)) 

-y J dr0 2 (r)(*(r,t)t*(r,*)p-p*(r,t)t^(r,iP9) 

Similarly to the case of one detector described in Sect. 
IIII1 we solve the master equation eq. (|19p by performing 
the transformation p(t) = F 2 {t)p(t)Fl(t) and 

^(r,t)=F^(t)^(v,t)F 2 (t), (20) 
where the operator F 2 (t) is defined as 



F 2 (t) = Te~ ei/2j, rft'/rfr'O 1 (r')* t (r',t')*(r',t')-e 2 /2/ t [ it'/dr'O 2 (r')1' t (r',t')*(r',t')^ 

I 



(21) 



The evolution of the rotated density matrix is thus given where 
by 



Kr,t) =e x j dr'ni(r / )*(r,t)p*t( r>t ) 

+e 2 J dr'n 2 (r')^(r,t)p¥(r,t), (22) 

The equation can be solved using perturbation theory, 
where we get an expression as in eq. (j9]) that includes 
correlation terms between the two detectors. 

The conditional probability distribution of counting m 
particles in one detector and n particles in the other one 
thus reads 



(-1)™+™ d m+n n 
P{m > n) = mini d\?d\- 2 Q 



Ai=l,A 2 =l 



(23) 



Q(Ai,A 2 ) = Tr(pT : e - AlIl - Aal2 :), 



and 



(24) 



T i (r)=e i J^ dtJ<KU(T)&{T,t)9(T,t). (25) 



The modified field operator that includes the absorption 
at the two detectors is given by 



*(r,t)= / dr'(r\e-^ t - to ^ n( - H °- ie ^ ni+a ^\r')^(r',to). 

(26) 
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V. DETECTION OF AN EXPANDING 
BOSE- EINSTEIN CONDENSATE 

Let us now illustrate the space and time dependent 
counting process by considering the counting statistics 
of a freely expanding Bose Einstein condensate. For 
simplicity, we consider a one dimensional system with 
a point-like detector located at a distance zq from the 
condensate. The detection time is of the order of the 
system dynamics, such that we calculate the full time- 
and space dependent generating function with the inten- 
sity given by eq. (1251) . We consider a point-like detector 
placed at a distance zq from the center of the cloud. The 
detector is modeled by a delta-function S(z ~ z ), such 
that the intensity eq. (j2"5)l reads 



1 = 



dt'&{zo,t')*{zo,t'), (27) 



where the time evolution of the operators ^>(r,t) is de- 
scribed by eq. (fT3|) . For the detection of a ID BEC at 
a point-like detector, the time evolution of the single- 
particle wave function is given by 



z,t)= dz'G(z,z',t,t Q )cl>(z',t ), (28) 



where 

G(z, z', t, t ) = ( z \ e -«t-t°)MH -iHe6) |^ ( 2 g) 

The counting distribution is then obtained from eq. 
([TO)) , which for the case of a condensate with TV particles 
reads 



p(m) 



(Ve / dtft(z ,t)$(z ,t) 



-Ne f dt<jJ(z ,t)<i>(zo,t) 

i 

(30) 



In the following subsection, we exactly solve eq. 
where we approximate the initial wave function <j)(z, 0) 
by a Lorentzian function, 



(z,0) = Vfe~ r|z| . 



(31) 



In Subsect. IV Bl we calculate an approximate solution 
obtained by the Born approximation. In Subsect. IV Cl wc 
compare the exact solution to the approximate solution 
as well as to the quantum counting formalism. 



A. Exact solution 

We calculate the counting distribution p(m) for an ex- 
panding BEC, where the detector is located at some dis- 
tance zo from the center of the cloud. We follow the 
treatment in [26| to derive an exact solution for the prop- 
agator eq. (|29|) that describes the whole system evolu- 
tion including the absorption at the detector. The sys- 
tem Hamiltonian is composed of two parts: the free par- 
ticle Hamiltonian Hq and a part corresponding to the 
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FIG. 1: Density |</>(zo, i)| 2 with respect to time (Fig a) and 
normalized intensity I/N with respect to r (Fig. b) for zq — 
(blue squares), zq = 0.1/T (green circles), zq = 0.4/r (red 
diamonds) . As the distance zo between the center of the cloud 
and the detector increases, the wave function spreads and 
long detector opening times are required to achieve non-trivial 
intensities. Parameters used: zo — 0, e = 1, N = 100, F — 10 5 . 



detection process, which acts as an imaginary potential. 
The propagator G(z,z',t) that describes the time evo- 
lution of the wave-function including the absorption at 
the detector can be written in an iterative way using the 
Lippmann-Schwinger equation, which for a point-like de- 
tector Sl(z') = 8(z' — zq) reads 



G{z,z',t) = G (z,z',t) - 
dt'G (z, z , t - t')G(z ,z', t'). 



(32) 



In the Appendix, we solve the Lippmann-Schwinger 
equation and show that the full propagator that describes 
the back-action of the detector is given by 

G(z ,z',t) = G (z ,z't) 

ime , ,,, . , me h . . , 

+—M(z-z + z'-z ,- — ,-t), 33 
h n m 

where M(z,k,t) denotes the Moshinsky function [28| 



M(z,k,t) 



dz'- 



exp ikz + i 



, , A z ~ z ') 2 



2t 



y/2mt 

- exp (ikz - i-k 2 t] erfc(e-" /4 ^^X34) 
2 \ 2 J \/2t 



The modified wave function <f>(zo,t) can then be calcu- 
lated by standard integration techniques using eq. (j2"51) . 
and the counting distribution p(m) obtained by eq. (|30l) . 
The counting statistics are determined by the time inte- 
gral over the square of the wave function (j)( z 0i t). Fig. Q] 
a shows the square of the wave function with respect to 
time for different distances zq. As the distance increases, 
the wave function spreads in time. In order to obtain 
non-trivial results, long opening times are required. Fig. 
[T]b shows the intensity of particles at the detector with 
respect to the opening time r for detectors placed at dif- 
ferent distances from the detector. 
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B. Born Approximation 

In the previous section, we obtained an exact solution 
for the time evolution of the wavefunction by solving the 
Lippmann-Schwinger equation (|32[) . In this section, we 
use the Born approximation in order to derive an expres- 
sion for <f>(z,t) in terms of the known propagator Go- In 
the second order approximation, we obtain 

G(z, z', t) = G Q (z, z', t)~ J dt'G Q {z, zo, t-t')G (z , z' , t') 

(35) 

This implies that up to second order, the solution to eq. 
(|2g| is given by 

4>[z,t) = <p{z,t) 

-^J dt'j dz'Go(z,zo,t-t')G (z ,z',t')(t>(z',to). 
m) — 

U! dt '^ M{z °'- ir ' m ' /m)) (36) 

Eq. (j36) describes the evolution of the wave function, 
where the absorption at the detector is taken into ac- 
count up to second order. We get the higher order Born 
approximation by writing eq. (1361) in exponential form, 

4>(z,t) = VfM(zo,~iT,ht/m) x 

e -e/hf* it' y/ M(* , -*r, tit' /m)/M(z ,-»r,fit/m) 

In Sect. IV C[ we show that the Born approximation de- 
scribes the situation more accurately than the quantum 
counting formalism. However, the effect of the absorp- 
tion is under-estimated. 




t (S) " T (5) 



FIG. 2: a) |0(O,t)| 2 with respect tot. b) Normalized intensity 
X/N with respect to r. We compare the exact solution (red 
triangles), second order Born approximation (blue squares), 
exponential Born approximation (light blue diamonds) and 
the quantum counting formalism (green circles). Both the 
quantum counting formalism and the Born approximation 
reach intensities that exceed the total number of particles for 
long opening times r. In the exponential Born approximation 
the intensity is bounded for large r, however, the asymptotic 
value exceeds the bound given by the exact solution. Param- 
eters used: z =0,e = 1,N= 100, F = 10 5 . 
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FIG. 3: Counting distribution p(m) obtained from exact solu- 
tion (red triangles), Born approximation (blue squares) and 
the quantum counting formalism (green circles). For short 
detection times (Fig. a) the approximations agree reasonably 
well with the exact solution. For longer detection times (Fig. 
b), the approximations are no longer valid. Parameters used: 
2o = 0,e= 1,N= 100, T = 10 5 . 



C. Effect of the Absorption at the Detector 

Let us now analyze the effect of the back-action of the 
detector on the field. From eq. (|30|) it is clear that the 
important quantities to study are the square of the wave 
function, <jf(zo, t)(j>(zo, t), its time integral, as well as the 
full counting distribution. We discuss the limits in which 
the quantum counting formalism and the Born approx- 
imation give valid results, and study the limitations of 
the approximative solutions. 

In Fig. [2k, we plot the square of the wave function 
with respect to time, and compare the exact solution 
to the solutions obtained by the born approximation and 
the quantum counting formalism. The exact graph corre- 
sponding to the exact solution decays more rapidly, as the 
absorption at the detector is considered. The Born ap- 
proximation underestimates the decay of the wave func- 
tion and thus the absorption, however, it describes the 
behavior more accurately than the quantum counting for- 
malism, where absorption is not considered. 



The effect is seen more clearly when studying the in- 
tensity of the field at the detector, which is given by the 
time integral e J*7 dt<f>* (zo,t)<j)(z ,t) (Fig. [5b). For short 
detection times, the exact solution and the approximate 
solutions coincide. As the detection time increases, the 
intensity of particles is overestimated both for the Born 
approximation and for the quantum counting formalism. 
Note that for long detection times, the second order Born 
approximation diverges, whereas the exponential Born 
approximation reaches an asymptotic value. 

Finally, we compare the counting distributions ob- 
tained by the exact solution eq. ([55)1 to the solution ob- 
tained by the quantum counting formalism and the Born 
approximation. The effect of absorption is clearly visi- 
ble in the counting distribution, where the approximate 
solutions deviate increasingly from the exact solution as 
the measurement time increases (Fig. [3]). The counting 
distribution calculated with the full formalism including 
detection is clearly different from the approximated so- 
lution. 
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VI. SUMMARY AND CONCLUSIONS 

We have derived a formalism for describing the count- 
ing distribution of space- and time dependent fields tak- 
ing into account the back-action of the detector on the 
field. We have illustrated the importance of the effect of 
the back-action for the free expansion of a Bose Einstein 
condensate. An approximate solution using the Born ap- 
proximation describes the behavior of the system more 
accurately than the quantum counting formalism. How- 
ever, for typical detection times of expanding BECs, the 
effect of absorption is under-estimated significantly both 
by the quantum counting formalism, as well as the Born 
approximation. We thus showed that for certain experi- 
mentally relevant situations, the full time- and space de- 
pendent formalism has to be applied. 
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The Laplace transform of the free propagator G(zo, z', s) 
is given by 

The inverse Laplace transform of G{z, z' , s) can be per- 
formed by standard methods [271, such that the propa- 
gator is given by 

G(z ,z',t) = G (z ,z't) 

+— M (\z-zo \ + \z'-zo ,— -,-t , A.6 
n n m 

where M(z,k,t) is defined in eq. (|54")l . The function 
(f>(z ,t) is then calculated by eq. (|2"8"j) . 

Note that for a detector placed at zq = in the center 
of the cloud, the counting formula can also be obtained 
by directly solving the time dependent Schrodinger equa- 
tion, 

- h 2 d 2 ~ 

ih(f)(z, t) = -—^(f>(z, t) — iHeS(z)(f)(z, t) (A. 7) 

2m oz 

with the initial condition given by 

4>(z,0) = VTe-^i. 



Appendix: Derivation of the Propagator Including 
the Back-action of the Detector 



We solve the Lippmann-Schwinger equation for the 
propagator G(z,z',t) that describes the time evolution 
of the wave-function including the absorption at a point- 
like detector f2(z') = S(z' — zq), 

G{z,z',t) = G (z,z',t) - 
dt'G (z, z , t - t')G(z , z', t'), (A.l) 



ft Jo 

where the propagator Gq{z, z',t) for the free expansion 
is given by 



, , m i\z - z' 

G {z,z ,t) = \/r-7rrexp( 



r\2 



2niht 



2ht 



(A.2) 



We perform a Laplace transform of eq. (IA.1I) and use the 
convolution theorem to get 

G(z, z', s) = G {z, z', t) - ^G (z 7 zq, s)G(z Ql z', s)(A.3) 

n 



From eq. $27\ , we observe that we are only interested in 
the propagator at z — z , such that 

rti i \ Go(zq,Z ,t) 

G(z , z ,s) = 77 r- -. (A.4) 

1 - Ti G {zo,z ,s) 



In order to solve eq. (|A.7j) . we first express it in terms 
of its fouricr transform. From eq. ([27]) is is clear that 
for a detector placed at zq = 0, we are only interested 



in 0(0, t) 



J_ oo dk(f)(k, t). The fourier transformed 



equation is a differential equation with variable coeffi- 
cients and can be integrated by standard methods [29| . 
We get 



4>(o,t) = Mt) 



/27rft 



K{t-t')4>(0,t')dt', (A. 



where <f>o (t) 



l r= 

/ 2ttK J- 



= Jdke^T(j)(k,0) and K(t,t') 



dke~ 



We take the Laplace 



/27T/IJ-00 y/2(t-t') 

transform of eq. (IA.8I) and use the convolution theorem 
[H 29.2.8] to obtain 



0(0, s) 



4>o(s) 



k(s)' 



(A.9) 



The term <^o(s) is calculated by using the residue theo- 
rem. Using the method of partial fractions the expression 
for 0(0, s) is written in the form 



verse Laplace transform is given by 



such that the in- 



8 



(fir + ime) 



(A.10) 
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